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We report neutron scattering measurements of the structural correlations associated with the 
apparent relaxor transition in Ki_ 2 ;Li 2 ;Ta 03 for x = 0.02 (KLT(0.02)). This compound displays a 
broad and frequency-dependent peak in the dielectric permittivity, which is the accepted hallmark of 
all relaxors. However, no evidence of elastic diffuse scattering or any soft mode anomaly is observed 
in KLT(0.02) [J. Wen et al., Phys. Rev. B 78, 144202 (2008)], a situation that diverges from that in 
other relaxors such as PbMgi/ 3 Nb 2 / 303 . We resolve this dichotomy by showing that the structural 
correlations associated with the transition in KLT(0.02) are purely dynamic at all temperatures, 
having a timescale on the order of ~ THz. These fluctuations are overdamped, non-propagating, 
and spatially uncorrelated. Identical measurements made on pure KTaOs show that they are absent 
(within experimental error) in the undoped parent material. They exhibit a temperature dependence 
that correlates well with the dielectric response, which suggests that they are associated with local 
ferroelectric regions induced by the Li’^ doping. The ferroelectric transition that is induced by the 
introduction of Li"*' cations is therefore characterized by quasistatic fluctuations, which represents 
a stark contrast to the soft harmonic-mode-driven transition observed in conventional perovskite 
ferroelectrics like PbTiOs. The dynamic, glass-like, structural correlations in KLT(0.02) are much 
faster than those measured in random-field-based lead-based relaxors, which exhibit a frequency 
scale of order of ~ GHz and are comparatively better correlated spatially. Our results support 
the view that random fields give rise to the relaxor phenomena, and that the glass-like dynamics 
observed here characterize a nascent response. 


I. INTRODUCTION AND COMPARISON WITH 
LEAD-BASED RELAXORS 

Relaxor ferroelectrics are technologically important be¬ 
cause they exhibit enormous room-temperature piezo¬ 
electric coefficients and sizable dielectric constants with 
relatively little hysteresisJ^ These materials are typi¬ 
cally characterized by a dielectric permittivity that peaks 
broadly with temperature and displays a frequency- 
dependent amplitude and position spanning at least 14 
decades in frequency^ This behavior differs markedly 
from the sharp, frequency-independent dielectric re¬ 
sponse observed in conventional ferroelectrics such as 
PbTiOa^ri The most studied relaxor ferroelectrics are 
the lead-based compounds having the general formula 
PhBOs, for which the R-site is occupied by one of 
two or more heterovalent cations. Two such relaxors 
are Pb(Mgi/3Nb2/3)03 (PMN) and Pb(Zni/3Nb2/3)03 
(PZN), both of which share the PbTiOs cubic perovskite 
structure^ii^ Understanding the origin of the unusual di¬ 
electric properties of these materials is a topic of consid¬ 
erable current interest. To help clarify the origin of the 
“relaxor transition” in the lead-based compounds, and 
to guide the development of new materials with simi¬ 
lar dielectric properties, it is instructive to compare the 


physical properties of PMN and PZN to lead-free sys¬ 
tems that display a similar frequency-dependent dielec¬ 
tric response. To this end we present a neutron inelastic 
scattering study of the relaxor Ki_a,Lia;Ta03 (KLT) for 
X = 0.02 (KLT(0.02)) and compare our results with those 
from prior studies of PMN-based and PZN-based relaxors 
as well as the undoped (non-relaxor) parent compound 
KTa03. 

Neutron inelastic scattering methods have played a 
historically significant role in shaping our understand¬ 
ing of ferroelectric materials. Studies of the lattice dy¬ 
namics, and specifically of the soft, zone-center trans¬ 
verse optic (TO) mode, have shown that a direct link 
can be made between the neutron inelastic scattering 
cross section and the bulk dielectric responseJ^ Early 
neutron scattering studies of several lead-based relax¬ 
ors, — which were later corroborated with optical 
measurements/^ discovered a strongly-damped, soft, TO 
phonon branch that was associated with the formation 
of local polar nanoregions (PNR) (elastic diffuse scat¬ 
tering). The unusual damping of the soft TO phonon 
branch, which was termed the “waterfall” effect, was 
initially considered to be a signature of relaxor behav¬ 
ior. However subsequent studies found similarly damped 
TO modes in heavily-doped 40%Pb(Mgi/3Nb2/3)03- 
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60%PbTiO3 (PMN-60%PT)ii and Pb(Zri_^Ti ^)03 
Neither of these materials is a relaxor, and both show a 
sharp, well-defined ferroelectric transition similar to that 
in pure PbTi 03 . Therefore the waterfall effect, though 
somewhat unusual for perovskites, cannot be regarded as 
a dynamic signature of the relaxor phase. 

Neither neutron nor x-ray scattering studies have 
found any evidence of a truly long-range ordered, struc¬ 
tural phase transition in PMN, as would be manifested 
by the splitting of a nuclear Bragg peak. There is, how¬ 
ever, a striking lack of consensus regarding PZN. Several 
groups have reported a structural transition in PZN from 
cubic to rhombohedral symmetryJ^^— But other groups, 
using high-energy x-ray scattering techniques, have pre¬ 
sented compelling evidence that this transition is con¬ 
fined to the near-surface region or “skin” of the crystal, 
and that the actual shape of the PZN unit cell within 
the interior of the crystal remains cubic at all tempera¬ 
tures The conclusion of a skin effect in PZN remains 
controversial]However similar observations have also 
been reported in pure PMN^ as well as various composi¬ 
tions of PMN and PZN doped with PbTi 03 i^‘ 2 ^i^ Most 
recently a skin effect has also been observed in the lead- 
free perovskite relaxor Nai/ 2 Bii/ 2 Ti 03 (NBT)i^ In all 
cases, the skin effect is observed in systems that possess 
strong, quenched, random electric fields. 

Despite the anomalous soft mode behavior and the 
controversial structural properties, neutron and x-ray 
scattering experiments on PMN and PZN have revealed 
two unique features that do appear to be experimentally 
linked to the relaxor state. These features provide the 
motivation for our study of KLT(0.02). First, both PMN 
and PZN exhibit strong, temperature-dependent, elas¬ 
tic (i. e. static) diffuse scattering that is polar in nature 
and highly anisotropic. This has been demonstrated by 
x-ray scattering measurements in zero and non-zero ap¬ 
plied electric field ^>21 neutron scattering studies 
and studies of the low-energy transverse acoustic (TA) 
phononsThe diffuse scattering is not purely static 
at all temperatures and in fact displays a slow relax- 
ational character on a GHz timescale that was identified 
using neutron spin-echo22 and backscattering22ii£ tech- 
niques4ir— In particular, the neutron spin-echo experi¬ 
ment by Stock et al. revealed a dynamic component to 
the diffuse scattering in PMN that they related to the 
peak response in the dielectric permittivityi22 A sum¬ 
mary of how the diffuse scattering may be a general prop¬ 
erty common to all relaxors, and not just to the lead- 
based systems, is given in Ref. [2^ The second unique 
feature of relaxors concerns the presence of so-called 
“columns” of inelastic scattering, broadly distributed in 
energy, that are located at the M-point {Q = {^, ^,0)) 
and i?-point {Q = (^, ^)) zone boundaries4i A similar 

zone-boundary soft mode was observed at the M-point 
with x-ray inelastic scattering methods in the disordered 
perovskite PZT..-4A The goal of our study of KLT(0.02) 
is to determine how universal these features are to the 
relaxor phase. 
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FIG. 1. Neutron inelastic scattering spectra from KLT(0.02) 
and KTaOa measured on the IT thermal neutron spectrome¬ 
ter located at the LLB. An underdamped TA mode is shown 
in panels a) and c) for both samples. Energy broadened relax- 
ational dynamics are observed in KLT (panel 6)) where the in¬ 
tensity decreases monotonically with increasing energy trans¬ 
fer. By contrast, KTaOs displays no such dynamics (panel 
d)) at the same wave vector and temperature. The horizon¬ 
tal bars represent the instrumental elastic energy resolution 
FWHM (full-width at half-maximum). 


KTa 03 is a stoichiometric, parent compound of KLT 
for which the ferroelectric-active, low-energy TO mode 
softens monotonically with decreasing temperature^^^Ti^S 
at least for all temperatures above ~ IK as discussed 
in Ref. IsH . No long-range ordered structural phase tran¬ 
sition, characterized by new, sharp (resolution-limited) 
Bragg peaks, has ever been observed. For these reasons 
KTa 03 is classified as an incipient ferroelectric. It is 
known that replacing K+ with Li+ on the perovskite A- 
site stabilizes a ferroelectric ground state in KLT for suf¬ 
ficiently large x > Xc- Microscopically, the much smaller 
Li+ ion shifts away from the high symmetry position 
along one of six [100] directions thereby creating a lo¬ 
cal electric dipole, which has been studied using NMR 
techniquesi2^ It has also been shown that Li+ doping 
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actually hardens the soft zone-center TO mode, which 
suggests that the transition to the ferroelectric phase is 
order-disorder in character and not displacive .^^'^^ How¬ 
ever, consensus on the precise value of Xc is lacking. In 
1981 Prater et al. measured the Raman scattering from 
various KLT samples as a function of temperature and 
concluded that they were consistent with the presence 
of a tetragonal or orthorhombic ferroelectric phase for 
X > O.Oli^ In 1987 Kleemann et al. performed optical 
birefringence, refractive index, and dielectric permittiv¬ 
ity measurements on KLT and proposed a crossover from 
a dipole-glass to ferroelectric phase at Xc = 0.022i^ At 
the same time dielectric measurements have shown that 
the peak in the dielectric permittivity of KLT is broad in 
temperature and strongly dependent on frequency across 
a range of concentrations 0.01 < x < 0.035 and thus sim¬ 
ilar to the dielectric response of PMN.— By analogy 
with the lead-based systems, the elastic diffuse scatter¬ 
ing in more heavily Li''"-doped KLT samples has been in¬ 
terpreted as resulting from short-range ferroelectric cor¬ 
relations or polar nanoregionsi^ These results are sup¬ 
ported by earlier neutron studies of the TA phonons in 
KLT(0.035), which employed an external electric field to 
observe the splitting of the TA modes resulting from a 
cubic-to-tetragonal structural phase transition.— More 
recently a neutron scattering study by Wakimoto et al. 
examined the lattice dynamics of a single crystal of KLT 
with X = 0.05 and, just as in KTaOa, found a zone-center 
TO mode that softens continuously on cooling and no ev¬ 
idence of any long-range ordered structural distortion.— 
Wakimoto et al. also observed elastic diffuse scattering 
that increased monotonically on cooling below 200 K, but 
it was present only near some nuclear Bragg peaks such 
as (210) and absent near others such as (100) and (200), 
and the diffuse scattering contours were elongated along 
(100); this situation differs markedly from that observed 
in PMN and PZN. It is quite likely that some of the dis¬ 
crepancies between the aforementioned studies may be 
the result of uncertainties and/or gradients in the nomi¬ 
nal values of the Li’*' concentrations. 

While there have been extensive studies of KLT for x > 
Xc, there has been comparatively little effort for very low 
Li’*' concentrations for which a well-defined ferroelectric 
transition is absent. For x < Xc the dielectric constant 
displays a frequency-dependent and broad peak consis¬ 
tent with a relaxor phase, but no evidence of any elastic 
diffuse scattering nor any anomaly of the soft TO mode 
characteristic of this phase has been observed— Neu¬ 
tron experiments on KLT(0.02) have reported a strong 
enhancement of the nuclear Bragg peak intensity below 
a critical temperature that roughly matches the temper¬ 
ature at which the peak in the dielectric response occurs. 
The Bragg peak enhancement is possibly indicative of a 
release of extinction and is similar to the behavior ob¬ 
served in PZN— However no Bragg peaks were observed 
to split or distort, suggesting that the unit cell shape 
remains cubic at all temperatures. 

The dielectric permittivity and crystal structure of 


KLT(a; < Xc) are very similar to those of PMN and PZN. 
It is therefore puzzling that no other scattering signature 
of the relaxor phase, like those observed in lead-based 
systems, are seen in KLT(a: < Xc). This begs the ques¬ 
tion of whether or not the dynamic features observed 
in the lead-based systems are common to other relaxor 
systems not based on PbTiOa. In this paper, we re¬ 
port neutron inelastic scattering results that demonstrate 
the existence of a dehnitive dynamic signature of local¬ 
ized (i.e. short-range ordered) ferroelectric correlations 
in KLT(0.02). We observe strong, overdamped fluctua¬ 
tions that achieve maximum intensity at a temperature 
close to that at which the dielectric constant peaks and 
also where the nuclear Bragg peaks display a large en¬ 
hancement in intensity. These fluctuations occur on the 
THz timescale, which is much faster than those observed 
in the lead-based PMN and PZN, and are poorly corre¬ 
lated spatially. A summary of these results is presented in 
Fig. [11 which shows underdamped TA phonons for both 
KLT(0.02) and KTaOs (panels a) and c)) and a broad 
relaxational mode in KLT (panel b)) that is absent in 
KTaOs (panel d)). 


II. EXPERIMENTAL DETAILS 

We use neutron scattering methods to study the struc¬ 
tural dynamics of KLT and to connect the physics to 
the dielectric response. Neutrons are uniquely well-suited 
to probe condensed matter systems over a broad range 
of momentum (length) and energy (time) scales with 
good resolution in both (typically Ag ^ 0.01 and 
AE ^ 1 meV or 0.25 THz). This capability nicely com¬ 
plements dielectric measurements, which can only probe 
fluctuations near |Q| = 0 but are sensitive to dynamics 
on much slower time scales ranging from GHz to mHz. 
Therefore in order to understand how the fluctuations 
observed with dielectric measurements are spatially cor¬ 
related and how they influence harmonic phonons on the 
THz timescale, it is essential to use neutrons, which are 
able to probe fluctuations throughout the Brillouin zone. 

The KLT(0.02) single crystal used for this study is the 
same as that investigated in Ref. [gT]. It has dimensions 
0.5 X 1 X 2cm^ and a room temperature lattice constant 
of 3.992 A. The Li concentration was determined by es¬ 
timating the concentration in the melt and then verified 
from the position of the peak in the dielectric permittiv¬ 
ity with temperature. To test whether or not the scat¬ 
tering observed in KLT originated from Li'*' doping, we 
performed an independent series of measurements on a 
pure KTaOs crystal. This crystal is the same as that 
used in Ref. HI. 

Most of the neutron scattering measurements reported 
here were performed on the IT thermal-neutron triple¬ 
axis spectrometer located at the LLB-Saclay reactor. 
The incident and final (scattered) neutron energies were 
defined through the use of a doubly-focused PG(002) 
monochromator and analyzer, which were used with open 
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collimation. For measurements on KLT(0.02) the final 
neutron energy was fixed to Ef = 14.7 meV and a highly- 
oriented pyrolytic graphite (HOPG) filter was used to 
remove higher-order neutrons from the scattered beam. 
For elastic Bragg peak measurements an HOPG filter was 
also placed in the incident beam to further reduce the 
possibility of higher-order neutron contamination from 
the monochromator. Our measurements were conducted 
primarily at low energies (< 12meV) where the strong, 
overdamped signal is present. To obtain a reciprocal 
space map of the lattice dynamics, we performed mea¬ 
surements on the NEAT cold-neutron chopper instru¬ 
ment, located at the HZB research reactor, using a fixed 
initial neutron energy Ei = 14.1 meV and the detector 
coverage on the neutron energy gain side. A similar tech¬ 
nique was used to map out the soft columns of scattering 
located near the zone boundaries in PMN.— 

Measurements on pure KTaOa were made using the 
IT spectrometer, where Ei was fixed to 36meV us¬ 
ing a Gu(lll) monochromator. A fixed and large in¬ 
cident energy was chosen to rednce backgronnd result¬ 
ing from higher-order nentrons incident on the same 
area and hence improve sensitivity to any diffnse “qua¬ 
sistatic” scattering near E = Q. The final neutron 
energy was scanned nsing the (004) reflection from a 
PG crystal analyzer, and the intensity was corrected for 
the i&TiiOanalyzer) factor reqnircd to acconnt for the 
change in resolntioni^ The sample was cooled nsing a 
closed-cycle cryostat. Additional measurements of the 
TA phonon dispersion were performed using the BT9 
thermal-neutron triple-axis spectrometer, located at the 
NIST Genter for Neutron Research, using a PG(002) 
vertically-focused monochromator and a PG(002) ana¬ 
lyzer. The horizontal beam collimation sequence used on 
BT9 was 40-80'-S'-80'-open {S denotes the sample posi¬ 
tion). An HOPG filter was inserted in the scattered beam 
to remove higher-order neutron contamination from the 
monochromator. The final neutron energy was fixed to 
Ef = 14.7 meV and the energy transfer E = Ei — Ef was 
scanned by varying the incident neutron energy. We note 
that all data for both IT and BT9 were corrected for con¬ 
tamination of the incident beam monitor as described in 
detail in the appendix of Ref. [13 and in Ref. [H. To verify 
consistency between the different experimental configura¬ 
tions we cross-checked the low-energy TA phonons shown 
in Fig. 13] a) and b). 


III. ELASTIC RESPONSE AND STRUCTURAL 
TRANSITION 

The structural properties of KLT(0.02) have been re¬ 
ported in detail by Wen et al. in both zero and non-zero 
applied electric field)^ but for completeness we review 
the zero-field properties relevant to the inelastic stud¬ 
ies discussed below. In the experiments described in 
Ref. [fill and here, no splitting of the Bragg peaks or any 
temperature-dependent anomaly in the lattice constants 



FIG. 2. The change in Bragg peak intensity at Q=(2,0,0) 
as a function of temperature measured on IT (LLB, Saclay). 
The open and filled circles correspond to cooling and warming 
sequences, respectively. A large hysteresis is observed and 
the intensity displays a peak on warming. The warming and 
cooling rates were set to 2K/min. 

is seen. This is consistent with the fact that the Li+ 
concentration x of our sample is less than the critical 
concentration Xc (discussed earlier) required for a well- 
defined structural phase transition to occur. 

We observe a very large increase in the Bragg peak 
intensity on cooling. This is illustrated in Fig. |2l which 
shows the Bragg peak intensity at Q={2flfi) as a func¬ 
tion of temperature. A large hysteresis of order 10 K 
is seen, and on warming the intensity displays a dis¬ 
tinctive peak near 50 K, which is close to the tempera¬ 
ture where the broad and frequency-dependent peak oc¬ 
curs in the dielectric permittivity. The complex tem¬ 
perature dependence of the elastic scattering shown in 
Fig. 13] is not expected for either a first or second-order 
phase transition and is similar to the loss of extinction 
in the Bragg peaks reported in some lead-based relaxors 
like PZN^ As noted in previous studies on KLT with 
small Li+ concentrations (Refs. Ifill) . we do not observe 
any temperature-dependent elastic diffuse scattering that 
decorates the nuclear Bragg peaks as measured in Ref. 
[ 53 . The high-energy x-ray study of Wen et al. did find 
evidence of diffuse scattering along (100) directions, but 
these were shown to be a result of the large energy res¬ 
olution associated with the x-ray measurement, which 
integrated over low-energy phonons. Our measurements 
here have all been performed with a monochromator and 
analyzer, which provide an elastic {E = 0) energy reso- 




IV. LOW-ENERGY PHONON DISPERSION - 
RESPONSE OF THE TAi PHONON TO LI+ 
DOPING 
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FIG. 3. Panels a) and 6) show the low-energy TAi phonon 
peak positions for KLT(0.02) and pure KTaOs, respectively. 
A constant-Q map of the low-energy phonons measured on 
NEAT is illustrated in panel c) at 300 K that spans two Bril- 
louin zones (divided by the dashed line). The negative values 
of the energy transfer E indicate that the experiment was per¬ 
formed on the neutron energy gain side. The region of high- 
intensity (shown in red) at E = 0 results from elastic incoher¬ 
ent scattering; the high intensity regions near Q = ( — 1, 2,0) 
and (0,2,0) are due to the low-energy TA and TO phonons. 
Of particular interest for this study is the weaker region of 
intensity that exists between the well-defined acoustic modes 
at the two zone centers, which displays a weakly varying in¬ 
tensity ridge, d) shows a constant energy cut illustrating that 
there is no strong momentum dependence to this weak scat¬ 
tering throughout the Brillouin Zone. 


lution of ~ 1 meV full-width at half maximum (FWHM). 
Therefore, we can state with conhdence that KLT(0.02) 
does not undergo a bulk, long-range ordered, structural 
phase transition, but it does exhibit an anomaly that re¬ 
sults in a loss of extinction (higher intensity) in the Bragg 
peaks on cooling to low temperature. The origin and the 
nature of this is discussed later in the text and related to 
the dielectric response. 


In this section we review the low-frequency lattice dy¬ 
namical response of KLT(0.02) and KTaOa as a function 
of temperature. The low-energy TA and TO modes have 
been the topic of several studies for KLT(0.02) (Ref. HH) 
and KTaOs (Ref. SiQ. We note that a detailed neu¬ 
tron scattering study of the TA phonons in a more highly 
doped sample of KLT(0.035) found clear evidence of a 
structural distortion through the use of an electric field.— 
We performed surveys of the low-energy TA and TO 
phonons, and our results are in agreement with previ¬ 
ous descriptions. 

Fig. [3] a) and h) shows the dispersions of the TAi 
modes (transverse acoustic phonons that are polarized 
along [100] and propagating along [010]) for KLT(0.02) 
and KTaOs, respectively, at 15 K and 200 K. In agree¬ 
ment with the classic mode-coupling analysis of Axe et 
al.^ we hnd that the TAi mode energy decreases on 
cooling for both compounds, which reflects the softening 
of the TO mode. And while the energy scales of the TAi 
phonon branches in KLT(0.02) and KTaOs are similar, 
there are two important features to note here. First, the 
long-wavelength modes in KTaOs are noticeably lower in 
energy. This is consistent with the mode-coupling sce¬ 
nario and the fact that Li-doping stiffens the TO branch 
in KTaOs— Second, the relative change in the TAi 
mode energy with temperature in KLT(0.02) is slightly 
less than in KTaOs. This suggests that the TO-TA mode 
coupling is weaker in KLT(0.02), which is again consis¬ 
tent with the Li-induced stiffening of the TO branch. 
This situation bears an interesting similarity to the case 
of PMN versus PbTiOs, which was studied by Stock et 
al. who demonstrated that the TO-TA mode coupling 
in PMN is extremely weak and thus cannot be the cause 
of the waterfall effect— These phonon data therefore 
conhrm the elastic results of Fig. [5] in that both show 
that the Li+ dopants have a signihcant effect on the bulk 
properties of the KLT(0.02) crystal. 

To search for other effects of the Li+ dopant on the 
lattice dynamics that might be relevant to the elas¬ 
tic anomalies discussed in the previous section, we per¬ 
formed a reciprocal-space survey of the lattice dynam¬ 
ics of KLT(0.02) using time-of-flight neutron scattering 
methods. Fig.[3]c) illustrates an intensity map generated 
by a series of constant-Q slices that were measured us¬ 
ing the NEAT direct-geometry time-of-flight spectrome¬ 
ter (HZB, Berlin) on KLT(0.02). The experiment was 
performed using the neutron energy gain-side (energy 
transfer E = Ei — Ef < 0) at 300 K because the en¬ 
ergy resolution is narrower near the elastic {E = 0) posi¬ 
tion and progressively worsens at higher energy transfers; 
this is a configuration that mimics that of a triple-axis 
spectrometer and better matches the characteristics of 
the cross section being measured. The region of high 
intensity around the elastic position results from the rel- 





















6 



FIG. 4. Temperature dependence of the change in scattered 
intensity at the zone-boundary M-point Q = (0.5, 1.5, 0) at an 
energy transfer of 5 meV. The data were measured by tuning 
the spectrometer to the M point and taking the lowest tem¬ 
perature as the background. The dashed line represents the 
Bose factor at each temperature scaled to fit the temperature 
dependence below 125 K. 


atively large incoherent cross section of potassium. The 
lines of high intensity at Q=(0,2,0) and (1,2,0) repre¬ 
sent strong scattering from the low-energy TA and TO 
modes previously characterized. An especially interest¬ 
ing feature is the significant inelastic scattering seen near 
the zone-boundary X-point (5=(-0.5,2,0), which extends 
from A = 0 to at least ^ —4meV, close to the top of 
the acoustic TAi phonon branch. While the scattering is 
most obvious at the zone boundary, where the acoustic 
modes lie at a higher energy (Fig. [3]a) and bj), the scat¬ 
tering is not well defined in momentum as evidenced by 
the extent of the contours with respect to the size of the 
Brillouin zone. This inelastic scattering is broadly dis¬ 
tributed in H and, in this respect, does not strongly re¬ 
semble the so-called “columns” of scattering reported by 
Swainson et al. in PMN using similar techniquesJi The 
fact that the scattering is highly extended in momentum 
implies the presence of a highly localized object in real 
space. To explore this broad scattering in more detail, 
we performed additional measurements using a thermal- 
neutron triple-axis spectrometer which we discuss in the 
following section. 



FIG. 5. Constant-Q = (0.5,1.5,0) scans measured with the 
IT thermal-neutron spectrometer with Ef = 14.7 meV. Pan¬ 
els a)-c) shows scans at 40, 75, and 200 K respectively. The 
solid line is a fit to the relaxational lineshape described in the 
text. The dashed line is a temperature-independent Gaussian 
function used to describe the strong incoherent elastic scat¬ 
tering. The horizontal bar represents the experimental elastic 
energy resolution (full-width at half maximum). 


V. OVERDAMPED LOW-ENERGY 
FLUCTUATIONS 

Motivated by the unusual inelastic scattering shown 
in Fig. [3] and the prior discovery of columns of spectral 
weight (i. e. scattering that is broadly distributed in en¬ 
ergy, but comparatively narrower in momentum) in PMN 
located at the M- and i?-point zone boundaries^ we 
performed additional measurements of the dynamics in 
KLT(0.02) using thermal-neutron triple-axis techniques. 
Fig. |4] shows how the inelastic scattering intensity mea¬ 
sured at the zone boundary M-point Q = (0.5,1.5,0) 
at an energy transfer of 5 meV changes as a function of 
temperature. The low-temperature value is assumed to 
be the background. The dashed line represents the Bose 
factor, which gives the temperature dependence expected 
for a simple thermal excitation. The measured temper¬ 
ature dependence is clearly inconsistent with that of the 
Bose factor. Thus we investigated this scattering further 
in momentum, energy, and temperature. 

Fig. [5] displays a series of constant-Q scans measured 
at Q = (0.5,1.5,0) at 40, 75, and 200 K. With increasing 
temperature a broad “quasielastic” component, indicated 
here by the solid line centered at FI = 0, is observed 
in addition to the strong, narrow in energy, incoherent 
elastic component, which dominates the total scatter¬ 
ing near E = 0. This quasielastic scattering contrasts 
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sharply with that from a well-defined harmonic mode 
such as a phonon, for which the lineshape in energy would 
be characterized by an underdamped peak centered at a 
non-zero energy. The observed excitations are therefore 
overdamped and display an energy linewidth (which is 
inversely related to the excitation lifetime) that is consid¬ 
erably broader than the intrinsic resolution of the spec¬ 
trometer. 

To parameterize these broad and overdamped fluctua¬ 
tions, we have fit the constant-Q scans to the following 
modified Lorentzian function, which is characteristic of a 
relaxational lineshape with a single energy scale: 

^(^) = Xo[’i(S) + 1] X-7^. (1) 

' + (f) 

For completeness we have also performed the analysis us¬ 
ing a damped harmonic oscillator, and this is presented in 
the appendix. We find consistent results between the two 
different lineshapes, and given the fact that the modified 
Lorentzian has fewer free parameters we use this line- 
shape for the remainder of the paper. Fitting was done 
using two parameters plus a temperature-independent 
Gaussian function centered at the elastic E = 0 position 
to account for the incoherent elastic scattering from the 
sample and mount. The parameter xo is an overall am¬ 
plitude and can be related to the real part of the suscepti¬ 
bility via the Kramers-Kronig relation. The parameter 7 
describes the linewidth of the modified Lorentzian func¬ 
tion and is related to the relaxational time via 7 oc 1 /r, 
and [n{E) + 1] is the Bose factor. The same lineshape 
has been used to describe the quasielastic scattering ob¬ 
served in the relaxor PMN (Ref. [63) as well as in various 
spin-glass systems (Ref. . 

The quality of the fits to the data using Eqn. 0] are 
illustrated by a series of representative curves in Fig. [5] 
as well as in Fig.[3and Fig.| 8 l which are discussed below. 
The solid line is a fit to Eqn. 0] and the dashed line rep¬ 
resents the effective background given by the Gaussian 
function used to describe the elastic incoherent scatter¬ 
ing. We see that the model cross section does an excellent 
job describing the data over a broad temperature range. 

This analysis employs only two temperature- 
dependent parameters, which are represented in 
Eqn. 0] by xo and 7 . The temperature dependence of 
both parameters is plotted in Fig. [ 6 ] The vertical dashed 
line in Fig. 0] is the average peak position of the real 
part of the dielectric permittivity (e') over the frequency 
range 100 Hz-1 MHz shown in Fig. 1 in Ref. 113 (which 
examined the same KLT(0.02) sample studied here). A 
clear peak is observed in xo at this temperature that is 
accompanied by a minimum of the energy linewidth 7 
at ~ 2meV. The peak in the susceptibility is also near 
where a large increase is observed in the Bragg peak 
intensity with decreasing temperature (Fig. [2]). We 
note that 7 measures the half-width of the scattering; 
therefore the quasielastic scattering always remains 



FIG. 6 . The parameter xo and relaxation rate 7 are plotted in 
panels a) and b) respectively. The parameters were extracted 
from constant-Q scans measured using the IT thermal triple¬ 
axis spectrometer at Q = (0,1.5, 0.5). The vertical dashed 
line indicates the peak position of the frequency-dependent 
dielectric response at xxx 200 kHz. The solid lines are fits to a 
Curie form (a) and a Arrhenius law ( 6 ). 


considerably broader than the energy resolution of the 
spectrometer, which is F = 0.55 meV (half-width at 
half-maximum). The response is not critical in nature as 
is characteristic of a long-range ordered structural phase 
transition. If this were the case, then a divergence of xo 
and 1/7 would be manifest (see, for example. Fig. 10 in 
Ref. IBOI . where a similar analysis was performed for the 
magnetic transition in Ba 3 NbFe 3 Si 20 i 4 ). The analysis 
shows that the broad quasielastic scattering shown in 
Fig. [5] is associated with dielectric properties and the 
broad peak in the dielectric response. To better confirm 
this point, we have performed measurements on KTaOs 
and find this scattering to be absent in the undoped, 
parent material. The anomaly of the energy linewidth 
near 70 K in Fig. [ 6 ] b) is reminiscent of the dynamics 
in the magnetic dynamics in the disordered multiferroic 
PbFei/2l'^bi/203'^ 

While our measurements at non-zero momentum 
transfer are, in principle, different from the dielectric re¬ 
sponse associated with |( 5 |= 0 , the connection between 
the peak in the dielectric constant and the anomalies ob¬ 
served XV. 60-70 K in 7 and xo discussed above suggest a 
common origin. We note, however, that the quasistatic 
scattering we observe is present across the Brillouin zone 
include the zone center and boundary. We cannot rule 
out the presence of zone boundary distortions and in the 
lead based relaxors, zone boundary and center phonon 
anomalies are observed at the same temperature4^ 

To compare our results directly we have ht (repre¬ 
sented by the solid lines in Fig. 0]) the high temperature 
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FIG. 7. Constant-Q scans performed at two wave vectors near 
the M-point Brillouin zone boundary located at (5=(0.5, 1.5, 
0). The values of the linewidth, 7 , in panels a) and b) are are 
3.0 ± 0.2 meV and 2.4 ± 0.2 meV, respectively. 
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dependence of xo to a Curie form xo cx 1/(T — ©vy) and 
r = tq exp (I7 //cb(T — ©vy)) (noting 7 oc 1/r and U is 
an activation energy) . The origins of such a “law” are 
described in Ref. We have included a “Weiss” tem¬ 
perature in both fits to xo and 7 to improve the fit and 
to account for the anomaly at ^ 60-70 K. We note that 
our data and temperature range are not able to distin¬ 
guish between a power law (r oc 1/(T — ©yy) suggested 
for spin glasses— and an Arrhenius form used here. 
We have chosen the exponential Arrhenius form to com¬ 
pare the energy scale of the fluctuations with dielectric 
measurements. Based on this analysis we derive a Weiss 
temperature of ©yy=38 ± 3 K, {7=187 ± 20 K, to=0.10 
± 0.02 ps. These values are very different from those de¬ 
rived from dielectric measurements where the energy bar¬ 
rier for the relaxational “site-hopping” process are U ^ 
1000-2000 K and tq 0.01 ps and associated with dy¬ 
namics of the Li'*' ions— The Curie temperature ©cvy is 
similar to the glass temperature measured with NMR and 
close to where a quadrupole split spectrum appears— 
However, given such disparate energy scales derived from 
Arrhenius fits, we conclude that the observed dynam¬ 
ics here are not directly associated with Li+ relaxational 
processes. To address this point and the origin of the 
dynamics, we now discuss the momentum dependence. 

The quasielastic scattering represented by a relax¬ 
ational lineshape is not only broad temporally (in en¬ 
ergy), but it is also spatially uncorrelated and hence dis¬ 
plays very weak momentum dependence. This is illus¬ 
trated in Fig.[3]by the broad scattering observed through¬ 
out the Brillouin zone. We explore this point further in 
Fig.H which shows a series of constant-Q scans near the 
(5=(0.5,1.5,0) zone boundary studied in detail above in 
Figs.[5]and[ni Scans for Q = {1 — q, H-g, 0) over the range 


FIG. 8 . a) and b) show constant-Q scans measured at 100 K 
for Q = (0.5, 1.5,0) and (1.5,1.5,0). The fitted linewidths are 
7 = 3.1 ± 0.2 meV and 4.5 ± 0.3 meV respectively, c) shows 
a plot of the half-width ( 7 ) as a function of Q^. The dashed 
line is a ht to Tick’s law for diffusion, and the solid line is a 
fit to a heuristic model for oscillating motion. 


oi q = [0.5,0.25] are shown at T = 100 K. The quasielas¬ 
tic scattering is clearly observed over this broad range in 
momentum transfer across a single Brillouin zone. The 
linewidth 7 also appears to decrease away from the zone 
boundary and at small values of [Qj (the values of 7 are 
given in the figure caption). The overdamped nature of 
this quasistatic scattering and the broadening with mo¬ 
mentum transfer does not represent a dispersion (as ob¬ 
served for the acoustic harmonic phonons as in Fig. ED. 

We next measured the scattering at different A/-points, 
shown in Fig. [ 8 ] a) and b), as well as at different \Q\, 
shown in Fig. [5] c). Fig. [3] shows the scattering near 
the Al-point Q = (0.5, 2, 0), whereas the temperature de¬ 
pendent measurements were made around the M-point 
Q = (0.5,1.5,0), shown in Fig. [5] Fig. | 8 ] shows analo¬ 
gous scans of the relaxational quasielastic scattering at 
T = 100 K. Panel a) and h) show scans performed at the 
M-points Q = (0.5,1.5,0) and (1.5,1.5,0), respectively. 
The different values of the fitted energy linewidth 7 are 
given with error bars in the figure caption. The quasielas¬ 
tic scattering is always present and broadens in energy 
with increasing |Q|. Measurements at equivalent Bril¬ 
louin zone positions and at different |Q| confirmed this 
trend. This proves that the scattering does not follow the 
symmetry constraints of the parent lattice and indicates 
the origin as highly localized objects. A discussion of 
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the relation between localized polar regions and the dif¬ 
fuse scattering and dielectric measurements is presented 

in Ref. [zl. 

The above analysis yields three key results. First, it 
demonstrates the presence of broad, quasielastic scatter¬ 
ing that peaks (measured through xo and 7 in Fig. [5]) at 
roughly the same temperature as does the dielectric per¬ 
mittivity and also where a sharp increase in the Bragg 
peak scattering intensity is observed. The quasielastic 
scattering is overdamped, and cannot be characterized 
as a well-defined dispersing harmonic mode. Second, 
the linewidth of the scattering increases with \Q\, but it 
does not follow the symmetry operations expected given 
the Brillouin zone fixed by the parent KTaOs lattice. 
Third, the activation energy (U) and the characteristic 
timescale (tq) are both an order of magnitude smaller 
than that measured in the dielectric response. Based on 
these three results, we suggest that the scattering is as¬ 
sociated with dilute, highly local droplets of dipoles, or 
ferroelectric regions, created by the introduction of Li'*' 
dopants, which enter randomly into the lattice. Given 
the ferroelectric nature of the localized object, we would 
expect a response at the zone center, consistent with di¬ 
electric measurements. Unfortunately the non-zero en¬ 
ergy and momentum resolution of the neutron spectrom¬ 
eter make measurements at the zone center problematic 
due to contamination from nearby acoustic modes. How¬ 
ever, by extrapolating the results obtained at non-zero 
wave vectors, we are able to infer that a response exists 
at the zone center that is consistent with the ferroelec¬ 
tric nature postulated here. The dynamics cannot be 
directly related to the hopping of Li+ cations given the 
large discrepancy between the activation energies derived 
here and from dielectric measurements. This is further 
confirmed by the fact that the integrated intensity of the 
relaxational “quasi static” component (Fig. [Ij is compa¬ 
rable to that in the low-energy acoustic phonons, which 
makes it unlikely to originate from the local hopping of 
a 2% dopant. These local dipoles are highly dilute and 
spatially uncorrelated as is evident from the absence of 
any well-defined, long-range-ordered, ferroelectric struc¬ 
tural transition. The comparison of the dielectric and 
neutron data performed on different timescales indicates 
a hierarchy of degrees of freedom in KLT(0.02) with dif¬ 
ferent energy scales. Such a scenario has been discussed 
in relation to the dynamics of glasses 

The scattering is unlikely to originate from a strictly 
zone-boundary distortion (such as tilting of the octahe- 
dra) because the temperature dependence of the neutron 
scattering linewidth and susceptibility display anomalies 
at the same temperature as do the \Q\ =0 dielectric 
measurements (Ref. [zl) and the Bragg peaks presented 
in Fig. [21 However, as noted by Swainson et al.^ the 
lead-based relaxors display both zone center and zone 
boundary phenomena. It should also be noted that the 
response observed here is proximate to zone-center diffuse 
scattering for slightly larger Li+ concentrations (Ref. IbOh . 
Besides the connection to the dielectric data, this behav¬ 


ior is also substantiated by calculations that sug gest that 
the distortion is a zone-center response (Ref. [T^ and also 
ultrasound measurements which probes the acoustic ve¬ 
locity in the limit as q approaches the zone center (Ref. 

US). 

We now discuss the momentum dependence of the en¬ 
ergy linewidth plotted in Fig. [5] c), which shows the half¬ 
width 7 as a function of |Qp. We emphasize that we ob¬ 
serve the fluctuations to be overdamped (the linewidth 
being larger than the energy position) at all momentum 
transfers and temperatures studied. The half-width 7 is 
well-described by a law for small momentum trans¬ 
fers, but this behavior seems to breakdown at larger mo¬ 
mentum transfers taken in higher Brillouin zones. Fol¬ 
lowing the study of localized dynamics in molecular sys¬ 
tems (Ref. HU), we fit our data to the following heuristic 
model, which has been applied to diffusing molecules for 
which the local oscillation time is much longer than the 
time for translational diffusion. 

Here we assume that the Debye term ~ 1 because 
the quality of the data in Fig. |5| does not support in¬ 
troducing a third parameter into the data analysis. The 
parameter 70 is inversely proportional to the oscillation 
time and a = Dtq, where D is the diffusion constant. In 
the limit of small momentum transfer Eqn. [2] reduces to 
Tick’s law (shown by the dashed line in Fig. | 8 ]c)), which 
states that the linewidth broadens in proportion to \Q\^. 
While the large momentum transfer data are somewhat 
scattered in Fig. |5|c), the data do show a broadening of 
the lineshape consistent with localized, and spatially un¬ 
correlated dipoles. The localized nature prevents them 
from having the same periodicity as the lattice and there¬ 
fore prevents the dynamics from being concentrated near 
the zone center. We speculate that this “diffusion” is 
tied to the fluctuations of locally ferroelectric distorted 
regions in a paraelectric background instead of the ac¬ 
tual physical hopping of Li+ cations. We emphasize that, 
given the overdamped nature of the fluctuations, this is 
not a dispersion as is defined for harmonic fluctuations 
(as shown for the acoustic phonons in Fig. ID. 

VI. CONCLUSIONS AND DISCUSSION 

We have presented a neutron inelastic scattering study 
of the low-energy lattice dynamics in KLT(0.02). This 
compound undergoes no long-range ordered structural 
distortion; instead it displays a strongly frequency- 
dependent dielectric permittivity that peaks broadly in 
temperature, which is very similar to that observed in the 
lead-based relaxor ferroelectrics. The main finding of our 
study is the existence of a broad quasielastic component 
that peaks in temperature near to where an anomaly is 
observed in the dielectric response. This scattering is 
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well described by a relaxational lineshape that is charac¬ 
terized by a single energy scale 7 oc l/r. 

The relaxational dynamics we report here are tied 
to local ferroelectric correlations induced by the Li+ 
dopants. We base this conclusion on several experimen¬ 
tal facts. First, the scattering we observe in KLT(0.02) 
is absent in pure KTaOa, which displays no dielectric 
anomaly at any temperature. Second, the spectral weight 
associated with the relaxational dynamics in KLT(0.02), 
probed through xo in Eqn. |4j peaks at the same tem¬ 
perature near to which the nuclear Bragg peak intensity 
suddenly increases and where a peak (albeit frequency- 
dependent and broad) is observed in the dielectric re¬ 
sponse. Third, the quasielastic scattering does not follow 
the symmetry of the KTaOa lattice; this indicates that 
it originates from local objects created by the dilute Li+ 
dopants. 

The presence of a broad relaxational component to the 
scattering in KLT(0.02) is corroborated by Raman mea¬ 
surements, which probe dynamics near the zone center. 
Prater et al. reported evidence of new excitations around 
50cm“^ ~ 6meV that broaden considerably with tem¬ 
perature in KLT(0. 014 ) 1 ^ Our results show that these 
excitations are broad in momentum and thus consistent 
with local (uncorrelated) objects. For larger concentra¬ 
tions of Li+ where a ferroelectric ground state is present, 
the local dipoles become spatially correlated and result 
in a soft TO mode that hardens at low temperatures and 
a subsequent splitting of the TA and TO phonons 

It is interesting to compare our results on KLT(0.02) 
with those obtained from the lead-based relaxor ferro- 
electrics represented by PbMgi/ 3 Nb 2/303 (PMN). PMN 
displays very strong and temperature-dependent diffuse 
scattering that is somewhat spatially correlated. While 
the diffuse scattering in PMN is also broadly distributed 
in momentum, it is also highly structured and has been 
described in detail in a number of x-ray and neutron stud¬ 
ies. The scattering observed in KLT(0.02) is compara¬ 
tively uncorrelated in momentum and better described by 
local, uncorrelated objects. The dynamics of KLT(0.02) 
are also significantly faster than those associated with the 
diffuse scattering observed in PMN. A neutron spin-echo 
study performed on PMN by Stock et ah^) found exci¬ 
tations with a timescale of order ^GHz. In KLT(0.02) 
we observe dynamics on the order of 1 meV or ^THz. 
The activation energies derived from the temperature 
dependence for PMN and KLT(0.02) are U ^ HOOK 
and 187 ± 20 K, respectively. There are therefore two 
key differences between the quasielastic diffuse scatter¬ 
ing present in PMN and KLT(0.02): (1) The scattering 
in KLT(0.02) is spatially uncorrelated in comparison to 
the strongly momentum-dependent and anisotropic dif¬ 
fuse scattering observed in PMN, and (2) the relaxational 
dynamics in KLT(0.02) are roughly two to three orders 
of magnitude faster than those in PMN. 

The relaxational or “quasistatic” component observed 
here in KLT(0.02) near the elastic {E = 0) line con¬ 
trasts strongly with the harmonic soft mode anomaly 


observed in KLT at larger Li'*" concentrations. These 
two rather disparate responses can be reconciled through 
the interpretation of the relaxational dynamics as a “cen¬ 
tral” component (i. e. one centered around A = 0 ) like 
that hrst discovered in SrTiOs by Shapiro et al^ and 
discussed theoretically by Halperin et alM- in terms of 
slowly-relaxing defects. Similar ideas have also been ap¬ 
plied to doped cuprate magnets^ where, in the small 
doping limit, slowly-relaxing defects produce a central 
component that gives rise to relaxational dynamics like 
those described for KLT(0.02). The lack of a strong mo¬ 
mentum dependence indicates weak doping. For larger 
doping concentrations the defects couple to the harmonic 
mode to produce a well-defined phase transition more 
characteristic of what is observed in the lead-based re- 
laxors. Although we draw this comparison, it is impor¬ 
tant to note that typical central peaks in the literature 
have energy widths that are comparable to the instru¬ 
mental resolution and phonon linewidths. The linewidths 
observed for KLT(0.02) are significantly broader, indi¬ 
cating a weak doping limit. The relaxational lineshape 
has some strong similarities to the dynamics observed 
near the zone center in KNb 03 , however the fluctuations 
in KNbOs are slower, having energy widths less than 
^ 1 meV, and comparatively well-correlated in momen- 

tum,S6 

The structural fluctuations in KLT(0.02) soften on 
cooling, yet no long-range ordered phase transition is 
observed with diffraction and there is no evidence of a 
soft TO phonon anomaly like that in PMN or PbTi 03 . 
Instead we observe only a strong increase in the Bragg 
peak intensity. We suggest that the low-energy, local, 
polar domains that we have associated with these fluc¬ 
tuations freeze near ^ 40 K, but that they do so in a 
random fashion given the lack of any strong momentum 
dependence. This likely creates a random local field that 
could then explain the history dependence and slow fre¬ 
quency dependence observed near the transition temper¬ 
ature. The Ising nature of the Li'*' displacement preserves 
the broad peak in the dielectric response and prevents 
the random field from entirely washing out the transi- 
tion^iiS^ This scenario is consistent with studies of model 
Ising magnets in a random field, for which strong hys¬ 
teresis and memory effects are also observed42i2£ This 
crossover from glassy (low-energy dynamics) to random- 
field behavior in KLT(0.02), which is characterized by 
a broadened dielectric peak and hysteresis, shares strong 
similarities with dilute ferromagnets including the lack of 
sharp propagating modes, the momentum dependence for 
7 , also the temperature dependencei^ Given that ran¬ 
dom helds play a key role in relaxor ferroelectrics, and 
given the close analogy with magnetic systems, we refer 
to KLT(0.02) as an incipient relaxor with properties that 
are more consistent with those of fluctuating glasses than 
with random fields. 

In summary, we have measured the dynamics associ¬ 
ated with local dipoles created by the substitution of Li+ 
for K+ in the incipient ferroelectric perovskite KTa 03 . 
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The low Li+ content of 2% is insufficient to disturb the 
cubic ground state structure. The dipoles are found to 
fluctuate on a timescale of ~THz, which is several or¬ 
ders of magnitude faster than that observed in the lead- 
based relaxors. The dipole dynamics, characterized using 
quasielastic neutron scattering methods, are broadly dis¬ 
tributed in momentum and thus consistent with uncorre¬ 
lated, spatially-localized objects. Given (1) the absence 
of any static diffuse scattering, ( 2 ) the local nature of 
the dipoles, and (3) the significantly faster dynamics, we 
suggest that KLT at these low Li concentrations not be 
classified as a relaxor, but rather as an incipient relaxor 
by analogy with the lack of ferroelectricity in KTaOa. 
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FIG. 9. Constant-Q = (0.5,1.5,0) scan measured at 100 K 
fit to a) a modified Lorentzian (Eqn. 3)) and b) a damped 
harmonic oscillator (DHO - Eqn. [5)l lineshape. The horizon¬ 
tal bar in panel b) represents the instrumental elastic energy 
resolution (FWHM). 


VIII. APPENDIX - DAMPED HARMONIC 
OSCILLATOR DESCRIPTION OF THE 
INELASTIC RESPONSE 

The low-energy scattering described in the main text is 
heavily damped and characterized by a linewidth that is 
always considerably larger than the instrumental energy 
resolution. This is evident from Fig. 6 where the full- 
width at half-maximum (2 x 7 ) of the scattering intensity 
is greater than 4meV. By comparison the instrumental 
energy resolution at the elastic line is measured and cal¬ 
culated to be 1.5 meV. Although the inelastic response is 
broad in energy, it does have a tendency to be peaked at 
very low energies (see Figs. 7 and 8 ). This situation is 
similar to that encountered in previous studies of amor¬ 
phous magnets (Ref. where two model cross sections 
were found to give equally good fits and physical interpre¬ 
tations of the data: 1) the “modified Lorentzian” (used 
here) and 2) the damped harmonic oscillator. Thus, 
while we favor the single modified Lorentzian descrip¬ 
tion, we will discuss the alternate description here and 
compare the results of this analysis to those presented 
in the main text so that the reader will be able to draw 
his/her own conclusions. 

The neutron inelastic scattering response I{Q,E) is 
proportional to S{Q,E), where S{Q,E) is directly pro¬ 
portional to the imaginary part of the susceptibility (y") 
times the Bose factor [n{E) + 1], 


I{Q,E)o^S{Q,E) (3) 

S{Q,E)^[n{E) + l]x''iQ,E) 


In the main text we modeled the broad quasielastic scat¬ 
tering using a modified Lorentzian lineshape that is char¬ 
acterized by a single energy scale defined by the linewidth 

7: 


x"{E) = xo -(4) 

l+(f) 

Here xo is related to the real part of the susceptibility. 
Alternatively, one can model the scattering in terms of a 
damped harmonic oscillator (DHO). This approach adds 
an extra parameter to the fit: the energy position 
The DHO lineshape has the following form: 


x"iE) = Xo 



^ E—hQp j ^ 


1 + 




. (5) 


Here we denote the half-width by T to distinguish it from 
that used in the modified Lorentzian lineshape in Eqn. 21 
A key requirement of all neutron cross sections is that 
x" be an odd function of energy; a property that both 
Eqn. in and [5] possess. In this Appendix we assess the 
results of the DHO analysis based on Eqn. [5] in terms of 
momentum transfer and temperature. 

Fig.inishows data from a constant-Q = (0.5,1.5,0) scan 
measured at 100 K. These data are fit to the modified 
Lorentzian lineshape (Eqn.2]) in panel a) and to the DHO 
lineshape (Eqn. [5]) ni panel b). It can be seen that while 
Eqn.mgives a good description of the data, Eqn.[5]fits the 
data slightly better because of the extra free parameter 
provided by the energy position hfto. 
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FIG. 10. The momentum dependence of the fluctuations near 
the M-point Q = (0.5,1.5,0) fit with the damped harmonic 
oscillator (Eqn. [Sjl lineshape described in the appendix. Pan¬ 
els a — d) show constant-Q scans at various momentum trans¬ 
fers. Panel e) illustrates the momentum dependence of the 
fitted energy position and /) shows the linewidth. The re¬ 
sults show that P ~ hflo, with both parameters exhibiting 
little momentum dependence throughout the Brillouin zone. 


We now investigate the momentum dependence of 
these parameters. Fig. [TOl illustrates a series of constant- 
Q scans measured near the M-point Q = (0.5,1.5,0). 
Panels a — d) show constant-Q scans measured along [110] 
that were fit to a DHO lineshape convolved with the spec¬ 
trometer resolution function (linewidth given by the hor¬ 
izontal solid bar shown in panel c). (We note that the 
anomalously low intensity point in panel d) near ~ 1 meV 
represents a change that is much sharper than the in¬ 
strumental energy resolution; it is also not reproduced at 
neighboring positions in momentum. Thus we attach no 
physical significance to this isolated data point.) Pan¬ 
els e) and /) display the momentum dependence of the 
energy position (Ttflo) 3 'iid linewidth (F), respectively. 
Within our experimental uncertainties, both parameters 



FIG. 11. A comparison of the temperature dependences of 
the linewidths F of the DHO lineshape (Eqn. [5]) and 7 of the 
modified Lorentzian lineshape (Eqn. m . 


are independent of q and furthermore Ti^l ^ F. The av¬ 
erage energy position is indicated by the vertical dashed 
line shown in panels a — d). From this we conclude that, 
independent of model (modified Lorentzian or DHO), 
the excitations probed here are overdamped and do not 
display the characteristics of propagating modes, which 
would manifest themselves as sharp, underdamped peaks 
in energy (like those shown in Fig. 1) that disperse with 
momentum. 

Finally, we compare the temperature dependence of F, 
derived from Eqn. [S] to that for 7 , derived from Eqn. |4j 
from which an activation energy was extracted in the 
main text. Fig. [TT] displays both 7 and F as functions 
of temperature. It can be seen that both parameters 
agree within error. (The DHO fit yields larger error bars 
because of the additional free parameter.) The value of 
the activation energy, which was discussed in the main 
text, is therefore independent of analysis. 

In summary, in this Appendix we have shown that 
the quasielastic fluctuations we observe in KLT(0.02) 
do not correspond to dispersive, well-defined, and un¬ 
derdamped excitations. We have also shown that the 
modified Lorentzian and DHO models yield a consistent 
answer for the activation energy. Given the extra param¬ 
eter of the DHO model and the ambiguity in interpreting 
the energy position parameter ?iDo due to its lack of q 
within the Brillouin zone, we favor the description pro¬ 
vided by the modified Lorentzian, which is characterized 
by a single energy scale. This is the analysis presented in 
the main text of the paper. We note that more complex 
lineshapes with multiple energy scales might also possi¬ 
bly describe our observations. An example is described 
in the context of spin fluctuations in superconducting 
YBaaCuaOe.ssa^*^^ or SrTiOg-^ But given the lack of 
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an underdamped peak in our neutron data, we have not found this analysis to produce underdetermined param¬ 
eters. 
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